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Abstract. - We present exact results obtained from Master Equations for the probability function 
P(y,T) of sums y — ^Yln=i^^ positions xt of a discrete random walker restricted to the set 

of integers between —L and L. We study the asymptotic properties for large values of L and T. 
For a set of position dependent transition probabilities the functional form of T) is with very 
high precision represented by ^-Caussians when T assumes a certain value T* oc L^. The domain 
of y values for which the ^-Caussian apply diverges with L. The fit to a ^-Caussian remains of 
very high quality even when the exponent a of the transition probability g{x) — \x/L\^ -\-p with 
< p << 1 is different from 1, although weak, but essential, deviation from the g^-Caussian does 
occur for a 7^ 1. To assess the role of correlations we compare the T dependence of P(y,T) for the 
restricted random walker case with the equivalent dependence for a sum y of uncorrelated variables 
X each distributed according to l/g{x). 



Introduction. — The central limit theorem states that appropriately scaled sums of 
independent random variables will be distributed according to a Gaussian [Tl|2]. The ran- 
dom walker is the prototype example of a stochastic Gaussian process [3l|4]. The standard 
random walker is characterized by transition constant probabilities, which are independent 
of position and time. Here we point out that for a certain class of position dependent transi- 
tion probabilities correlations arise, which leads to deviations away from Gaussian behavior. 
There exist already a large amount of evidence, which points to g'-Gaussians as the relevant 
high quality approximates for the functional form for the distribution function in a range 
of cases where correlations play an essential role [5 . The evidence for the relevance of the 
Gaussian is however often derived from numerical experiments in which fluctuations lim- 
its the accuracy and therefore the precision of the fit to the (^-Gaussian form (analytical 
exceptions to this frequent difficulty can be found in [6 ). Moreover we believe the random 
walk example we discuss here to be highly generic. It is related to e.g. particles moving in 
a confining potential or to more branching processes subject to resource limitations. 

Here we present an investigation of the sum y = Ylt=i of positions Xt passed through 
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by a Restricted Random Walker (RRW). The underlying stochastic process is sufficiently 
simple to allow exact numerical solution of the Master Equation (ME) for the probability 
distribution P{y^T). This ensures a very high precision fit to the g-Gaussian form and 
thereby a very accurate determination of the relevant parameters. We find that a broad 
range of transition probabilities for the random walker leads to g'-Gaussians with q param- 
eters depending on transition probabilities. Since the ME can be easily handled in exact 
numerically form, the RRW model is an excellent laboratory for understanding the condi- 
tions under which sums of correlated random variables are distributed as g-Gaussians: 

= I P(0) [1 - - q)y^] ^ for - q)y^ < 1 
[ otherwise 

where q < 3 and /3 > are parameters (for q > 3 normalizability is lost). As q ^ 1 the 
function P{y) approaches the Gaussian. 

Restricted Random Walk Model. — We consider a one dimensional symmetric 
random walker confined to the integers between —L and L. The motion of the walker is 
controled by the following time evolution 

{xt with probability g{x)/2 

Xt — I with probability g{x)/2 (2) 
Xt with probability 1 — g{x). 

We concentrate on the following form 

|-^| l|, (3) 

with reflective boundary conditions: If Xt-\-i > L {< —L) we let Xt+i ^ Xt+i — 1 (+1). 
We find numerically that the first return time (defined as the time elapsed until the walker, 
who leaves its x = position, returns to the zero position again, and note we do not 
include walkers who remain at x = for all times) distribution for these RRW behaves 
asymptotically like P{T) ^ ^ with r = 2, i.e. different from the exponent r = 3/2 for 
ordinary RW. We study the sum y = X]t=i in the limit p ^ for values of the exponent 
a = 0.75, 1 and 1.25. For p ^ 1 and L ^ oo the process reduces to the ordinary random 
walk. 

The highly restrictive nature of the RRW is clearly seen from the 6 trajectories shown 
in Fig (1] in the case L = 120 (with a = 1 and p = 5 • 10~^). For comparison we present 
the trajectories of an ordinary random walk on Xt-\-i = ± 1 with probability 1/2 and Xt 
confined to — L,— L + 1,..,L — 1,L. The figure shows that the vanishing transition probability 
g{x) near x = makes the RRW non-ergodic leaving most of the phase space empty. It is 
straight forward to derive a Master Equation for the distribution P(x, t) 

Px{x, t + 1) = Px{x, t) + ]^g{x - l)Px{x - 1, t) + ]^g{x + l)Px{x + 1, t) - g{x)Px{x, t) (4) 

subject to the appropriate boundary conditions at \X\ = L. The insert in Fig. [1] exhibits 
a solution of the equation. We will discuss P(x, t) in a more extended publication, here we 
now turn our attention to another distribution. 

Let P{y^ X, T) denote the probability that X^^^ Xt = y and xt = x. The time evolution 
of this simultaneous probability is controlled by the following ME 

P{y,x;t^l) = P{y,x;t)^ ^ 

AG{-1,0,1} 

[W{y, x;y-{x-A),x- A)P{y - (x - A), x - A; t)) 

-W{y + X + A, X + A; x)P{y, x; t)] . (5) 

p-2 



Restricted random walk model as a new testing ground for the applicability of g'-statistics 




Fig. 1: Six representative trajectories for L = 120. Main panel: The restricted RW model (p = 
5 • 10~^). Lower Inset: The standard RW model (p = 1). Non-ergodic behavior of the restricted 
RW model can easily be seen. Upper Inset: The time evolution of P(x, t) at three different t values 
(t = 1000, 5000 and 10000 from bottom to top). 



The transition probabilities W only depend on x and A. We have 
W{y^ x;y — {x — A), x — A) = w{x — A, A) 

and 

W{y + X + A, X + A; x) = w{x^ A) 

where 

w(z A) - { if A = ±1 

By substituting g{x) we obtain the following simple equation 

Piy, x;T+l) = ^[g{x + l)P{y -{x + l),x + 1; T) 
+gix - l)Piy - (x - l),x - 1;T)] 

+ {l-g{x))P{y-x,x;T) (7) 

The relevant boundary conditions are straightforward but lengthy to write down. 

We now investigate the functional shape of the distribution P{y, T) = P{y, x, T) for 
different values of a, L and T, for the typical value p = 5 • 10~^. In Fig. |2] we plot a typical 
case from where the perfect agreement between the exact and simulation results is evident. 
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Fig. 2: Exact and simulation results of the case a = 1, L = 120 and p = 5 • 10~^. It is clearly seen 
that the probability function P{y,T*) obtained from simulations is completely in accordance with 
the exact results. The number of experiments used in our simulations are 2x10^. 



Fig. [3] is concerned with the case a = 1 for different values of L. For fixed L we determine the 
value T* for which an optimal fit to a g'- Gaussian is possible. For T > T* values, each curve 
will start to exceed the g-Gaussian tails before the fast drop region, due to the finite-size, is 
being achieved. It is also interesting to analyze the relationship between T* and L. We find 
that T* L^, a behavior identical to the ordinary scaling that relates time and distance for 
diffusive processes. 

Since we have an exact numerical solution we can investigate with great accuracy the 
nature of the convergence to the g- Gaussian as we increase the domain L of the random 
walker. In Fig.|4]we demonstrate that as L is increased a trajectory in the T*-q-p parameter 
space exists along which P{y^T*^L) becomes increasingly well described by a g-Gaussian. 
The Figure contains the scaling combination Y = lng[P(?/, T)/P(0, T)]/[/3s^], where s = 
yP(0), \nq{x) = {x^~^ — 1)/(1 — q) is the g-logarithm and the scaling parameter Smax is 
defined as the s value of each L for which Y significantly starts to deviate from the — 1 
line, namely when + 1| > 0.004. If the dependence on y is exactly g-Gaussian, we would 
have F = — 1 for all y. An appropriate scaling of x-axis yields a clear data collapse. For 
all values of L we observe the deviation from —1 to be no more than a few parts in a 1000 
and as L increases the curves indeed approaches the line F = — 1 for large values of the 
argument y. The oscillations about the Y = —1 curve exhibit a subtle dependence on L. 
Careful inspection of the top panel in Fig. |4]reveals that for increasing values of L the curves 
actually approach the Y = —1 line for both small and large values of the argument s/smax- 
We therefore believe that asymptotically the distribution P{y^T*) indeed becomes very well 
described by the q- Gaussian functional form. It is unfortunately numerically impossible for 
us to reach very large L- values. 

This suggests the distribution P{y^ T, L) asymptotically is described by the g-Gaussian 
form if one let the pair (L, T) vary appropriately. We localize these very precise q and {3 
values so that the curves are as symmetric as possible along the —1 line. Using the q values 
given in Fig. HI we obtain an exponential dependence on L from where one can predict the 
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Fig. 3: The case a = 1 for L = 120, 140 and 160 with p = 5 • 10~^. The main panel shows the 
probability function P(y,T*). The center of the function is shown in detail in the left inset. The 
time T* is chosen to optimize the fit to the ^-Gaussian. The scaling of T* is given in the right inset. 



asymptotic value lim^^^oo q{L) = ^oo — 2.351, which is evident from Fig. [5l It is interesting 
to note that for values 5/3 < q < 3 the variance diverges. So in this respect the distribution 
behaves similarly to e.g. the Cauchy-Lorentz distribution, which corresponds to q = 2. 
Diverging variance is of course a common feature in complex systems of distributions with 
power law tails. 

Next we consider the effect of changing a to values different from one. Although it is 
still possible to tune T* so that the distribution P(?/,T*) is very close to a g-Gaussian, the 
high resolution Y plot given in Fig. [6] now shows that the order of deviation from straight 
horizontal line through Y = —1 grows significantly whenever a ^ 1. 

The role of correlations. — One might perhaps wonder to what extent the observed 
deviation from ordinary Gaussian behavior is caused by the peculiar shape of the probability 
distribution of the individual terms Xt in the sum y = XltLi ^t- To check this we solved the 
Master Equation for the probability distribution for y in the uncorrelated case where all the 
individual terms in the sum are drawn independently with probability 

= ^TTp' 

for X G {— L,— L + 1,..,0,..,L} and A/" the normalization factor. The motivation for this is 
simply that for the RRW considered above a term will appear in the sum y a number of times 
roughly given by l/g{x). In Fig. [71 we show that when the terms are uncorrelated the sum 
converges towards an ordinary Gaussian. We note that the uncorrelated distribution P(y, T) 
for small values of T does resemble a g'- Gaussian in the region of small y values. However, as 
T is increased the functional form rapidly changes towards the ordinary Gaussian in stark 
contrast to the correlated case (left panel in Fig. [7]) where the P{y^ T) grows towards the 
g-Gaussian as T is increased up to very large values of T. For the uncorrelated sum no 
trajectory (L^T^) which for L ^ oo takes one to the g'-Gaussian exists. 
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Fig. 4: The case a = 1 and p = 5 • 10~^ with L values between L = 90 to 160. The lower panel 
shows the data collapse when the x-axis is appropriately scaled. The upper panel shows a zoomed 
region around Y — —1 line. In the inset the scaling of Smax with L is given. The straight line is 
Smax AL^ with A 0.2726 and C 1.7828. 



Conclusions. — We have presented the hitherto most simple setting in which q- 
Gaussians control asymptotic behavior. We conclude that the Gaussian behavior is brought| 
about by the strong correlations and the high reluctance for the walker to move away from 
the central region of its domain. 

The numerical exact solution of Master Equations allows us to present high precision 
data for the probability function of sums of correlated random variables derived from a 
restricted random walk (RRW) with position dependent transition probabilities. When 
the range of the walker L and the number of terms in the sum T is scaled according to 
T = 1.54L^, g-Gaussians are observed over an increasingly broad interval. For non-linear 
transition probabilities we are able to identify a subtle oscillatory behavior away from the 
pure Gaussian form. Given the relative simplicity of the RRW it appears likely that the 
relation between transition probability and the value of q and the existence of oscillatory 
corrections to the Gaussian asymptote can be unraveled analytically. The RRW model 
presented in the present letter promises this way to significantly increase our understanding 
of the mechanisms responsible for the often encountered g-Gaussians. 

The very weak dependence of q on T'^{L) and the subtle oscillations in the a = 1 case and 
the more essential oscillations present for a 7^ 1 indicates that the true exact mathematical 
asymptote might not strictly be g-Gaussian but rather some functional form resembling a 
g'- Gaussian to a high degree of accuracy. It is only because we have numerically iterated the 
Master Equations exactly that we are able to identify this very slight difference. In studies 
relying on simulations or observational data the accuracy may not be sufficient to resolve 
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2.351- q = 0.015 exp[- (L/308.7) ] 



linear regression 
(r^=0.99912) 
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Fig. 5: Linear- log representation for the L dependence of q values. This exponential dependence 
suggests an asymptotic value around qoo — 2.351. 
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Fig. 6: Y plot of cases a = 0.75, a = 1 and a = 1.25 for a representative L value. Whenever a / 1, 
increasing order of deviation from —1 line is evident. 
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y P(0) y P(0) 

Fig. 7: Comparison of the T dependence for fixed values of L with p = 5 • 10~^ for the correlated 
(left panel) and the uncorrelated (right panel) P{y,T) distribution (see text for the details). 

these details and one would conclude that a Gaussian is an excellent approximation of the 
observed behavior. 

Let us recall that g-Gaussians have been found previously for more complex processes 
than the random walk to be able to provide very high quality approximations to relevant 
distributions. The case q < 1 was considered in [8l[9] where the authors found for two (scale- 
invariant) probabilistic models that the large-size limiting distributions are amazingly close 
to g-Gaussians, but are not exactly g'-Gaussians [10 . The work in Ref. [6 provides an 
analytic example of large-size limiting distributions that are g-Gaussians. We stress that 
even if g-Gaussians are not always the exact analytic form of the probability distributions 
in question, it is highly intriguing why they provide such exceptionally high accuracy ap- 
proximations in a large number of cases where correlations are sufficiently strong to make 
the central limit theorem inapplicable. 
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